Non-adiabatic charge pump: an exact solution 
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We derived a general and exact expression of current 
for quantum parametric charge pumps in the non-adiabatic 
regime at finite pumping frequency and finite driving am- 
plitude. The non-perturbative theory predicts a remarkable 
plateau structure in the pumped current due to multi-photon 
assisted processes in a double-barrier quantum well pump in- 
volving only a single pumping potential. It also predicts a 
current reversal as the resonant level of the pump crosses the 
Fermi energy of the leads. 
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Parametric charge pump is a device that drives the flow 
of a DC electric current by a time dependent variation of 
two or more device parameters. The quantum paramet- 
ric pump has received consid^ahk^l^lir^hDth exper- 
imentally and theoreticaUvyMMMElWa. A very 
useful and clear transport theory on quantum parametric 
pump was provided by BrouwerEJ for the adiabatic regime, 
i.e. in the lu — > limit, where lu is the frequency of driv- 
ing force of the pump. Since then there have been several 
attempts devoted towards the understanding of paramet- 
ric pumping at finite Jfcepuency, i.e. for thqjjapjje general 
non-adiabatic regime&H'Eil These theoriescloEj were all 
based perturbative schemes either in terms of pumping 
frequency or in pumping amplitude. Recently, a Floquet 
scattering matrix theory was developed!! 2 ] and an exact 
solution obtained for an oscillating double-barrier pump 
structure. Despite the progress, so far a general and non- 
perturbative theory for quantum parametric pump has 
not been reported. Such a theory will be very important 
to establish a general and unambiguous physical picture 
of quantum charge pumps for the non-adiabatic regime, 
and provide new features not accessible by perturbative 
theories. It is the purpose of this Letter to present the 
exact solution for the current delivered by a quantum 
parametric pump at both finite driving frequency and am- 
plitude. 

Our theory is based pii Keldysh non-equilibrium 
Green's functions (NEGF)B and is non-perturbative. It 
reproduces not only all previous theories in their perspec- 
tive limits, but also reveals qualitative new and interest- 
ing transport properties of quantum parametric pump. 
Our theory is also applicable to the case where the driv- 
ing potential is in the interior of the pump. Because of 
the generality of our theory, it is now possible to inves- 
tigate situations beyond the established physics of adi- 
abatic pump. As an example, we investigate a peculiar 
case of a quantum-well charge pump driven by only a 



single parameter: such a pump is only possible in the 
non-adiabatic regime. In particular, if the pumping po- 
tential is located at the barriers of the quantum well, 
the pumped current exhibits remarkable plateaus due to 
a multi-photon assisted pumping process. The width of 
the plateaus is determined by Tilu while their height vary 
depending on the pumping amplitude and the coupling 
strength between the pump and the leads. By tuning 
a gate voltage, the resonant level of the quantum well 
crosses Fermi energy of the leads and a pumped current 
reversal is observed. On the other hand, if the pump- 
ing potential is at the center point of the quantum well, 
due to the spatial symmetry the pumped current van- 
ishes exactly. Moving away from this symmetric point, 
the pumped current increases rapidly by several orders 
of magnitude and similar plateaus in pumped current are 
observed. 

Now we derive the pumped current in the non- 
adiabatic regime. Assuming that the time-dependent 
driving potential has the following form H'(t x ) — 
(V/2)exp(iu;ti) + (V*/2)exp(~iujti) where V is the ef- 
fective pumping potential profile. Neglecting interaction 
between electrons in the ideal leads labeled by L and 2, 
we start from the expression of time averaged currently 
(h = 1) in standard NEGF theory, 
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q j T dtf dt x Tr[G r (jt,h)Jl<(tut) 
G<(t,t 1 )S a a (t u t)+c.c.} (I) 



where r = 2t:/lu is the period of pumping cycle. Using 
the following double time Fourier transform 

&{t 1 ,t 2 ) = I dE2 e- iElti+iE ^G'y(E 1 ,E 2 ) (2) 
J 2ir 2tt 

with 7 = r, a, <, Eq.(0) becomes, 
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Tr{[G r (E 1 ,E 2 )-G a (E 1 ,E 2 )} 



2Nt J 2tt 2tt 

E < (E 2 , E x ) + G< (E 1 , E 2 ) [£« (E 2 , E x ) - V a {E 2 , E x )] } 

-^^Tr[(G''-G a )£< + G< (££-££)] (3) 



where we have extended the integration range for dt in 
Eq.(@) to [-Nt, Nt] with N— > 00. Note that in Eq.(§) 
the matrix G 7 has element G 7 (E±, E 2 ). It is also easy 
to show that G r - G Q = -iGTG". Using the fact that 
G< = G''S<G a and ^(E 2 ,E 1 ) = 2tt5{E 1 - E 2 )YT {E x ) 
(no external bias), we obtain 
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T(E 2 )G a (E 2 ,E 1 )](f(E 2 )-f(E 1 )) (4) 



we obtain 



Using the identity 

= 2iVr / dE/(2n) 

I7i J 



(5) 



we finally arrive at 



r dF +oa 
7 ^ E Tr[r a (E)G r (E 7 E- 



n— — oo 



a (2iVr) 2 

T n (E)G a (E + ncu, E)] (f n (E) - f(E)) (6) 

where T n (E) = T(E + nu) and f n (E) = f{E + nuj). Now 
we proceed to calculate G r (E, E + nuo). From the Dyson 
equation, we have 

G r {t,t') =G 0r {t,t')+ ( G r (t,t 1 )H'(t 1 )G 0r {t 1 ,t')dt 1 



(7) 

where G 0r is the equilibrium Green's function. Taking a 
double-time Fourier transform of Eq. (^) , we find 

G r (E 1 ,E 2 ) = 2irG 0r (E 1 )5(E 1 ~ E 2 ) 

+ / ^G r (E 1 ,E 2 + E)H'(E)G Qr (E 2 ) . (8) 

Since H'(E) = tt[V5(E + w) + V*5(E - u)], we obtain 

(r{Ex,E2) = 2TtG 0r (E 1 )6(E 1 - E 2 ) 
+ [G r {E ll E 2 - ui)V + G r {E 1 ,E 2 + oj)V*}G 0t '(E 2 )/2 . (9) 

Using this equation, the general expression of G r (E, E + 
ncu) can be obtained. To simplify notation, we use the 
following abbreviations G r n {E) = G r (E,E + nui) and 
G° n r (E) ee G 0r (E + nu). We then have 

G r = 2TrG a r S{0) + {G\V* + G r _ 1 V)G° r /2 

G[ = {G r 2 V* + G r V)Gl r /2 

G r _t = (G r Q V* + G r _ 2 V)G°T 1 /2 (10) 

If we restrict ourselves for single photon process only, we 
can neglect G 2 in Eq. (|To|) and obtain 



G r (E) = 



2tt(5(0) 



(11) 



where 2tt5(0) = J dE = 2Nt. Here the self-energy 
H\{E) ee (yG\ r V* + V*G°T 1 V)/4. 

The exact solution of G r n by including all photon pro- 
cesses can be obtained by iterating the following equation 
obtained from Eq.(||), 



Gr fiT jr/^iOr 1 



(13) 



where a r n = [(a n , a„), (a n +i, a n+ x), ...] with a n = 
iVG% +1 /2 and d n = iV*G° n r /2. The continued fraction 
[(ai, ax), {a 2 , a 2 ), ...] is defined as 



[(ai,Oi), (a 2 ,a 2 ), ...] = l + Oi- 



1 a x (14) 



1 + a 2 7— — a 2 

From Eq.©, we obtain G\ = G r Q V Gf / (2a\) 
and G r _ x = G r V*G°i: 1 /{2f3L 1 ). Here fi r _ n = 
[(6_ nj 6_ n ), (6_n_i, 5_ n _i),...] with fe_„ = iV*G^ n _j2 
and 6_„ = iVG°J n /2. Substituting the expressions of 
G±! into Eq. (|l0|) , we have 



G r (E) 



2tt(5(0) 



[G°r(E)]-i-Zr(E) 



with 



S r = VG 0rJ_y* + V * G 0r _J_y 

1 Aa\ ~ l W-x 



(15) 



(16) 



Once Gq is obtained, G r n can be calculated recursively 
from Eq.(|lJ). Eqs.(§), @ and © form the 

central result of this paperEJ. Note that G° r is the 
Green's function for the process of absorbing (emit- 
ting) n-photons whereas G r n is the renormalized n-photon 
Green's function containing the contribution of multi- 
photon process. From Eq.(13), we see that G r n is at least 
of the order of V n . For this reason, as we will see below 
that in the resonant tunneling regime, the summation 
in Eq.(^J) converges rapidly even in the strong pumping 
regime. Using the fact that GTG" = G a rG r , it is 
easy to show ^2 a I a — 0, i.e., current conservation for 
parametric pumping. 

Eq.(^) reproduces all previously known results of para- 
metric pumping. We give two examples here. First, in 
Rcf. O, the pumped current is obtained using perturba- 
tion theory for finite pumpinffiJjrcqucncy up to the second 
order in pumping amplitude!!^!. This result is recovered 
by only keeping the terms n = ±1 in Eq.(H) to obtain 
G\ = G° r VG° 1 r /2. Hence Eq.(|) becomes, 

/J 171 
— ]T Tr [r a G° r VG° n r r n G a n a V*G a a ] (/„ - /) 
71 n=±l 

(17) 

which is exactly the same as that of Ref. |ll|, |lj. Sec- 
ond, we consider the well-studied adiabatic regime lo — > 
where the instantaneous approximation for the Green's 
function is appropriate. We transform Eq.(^J) to the 
Wigner representation using 



G r n = (G r n+1 V* + G r n ^V)G°:/2 



(12) 



G~i(tx,h) 



dE 
2^' 



-iE(tx-t 2 ) syy 



G 1 {E,T) 



(18) 
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where T = (t 1 +t 2 )/2 and 7 is, the NEGF in the Wigner 
representation. We then havetj 



G r n = J dTQ r (E, T)e 



■inuT 



(19) 



Substituting Eq.(|19|) into Eq. 
term, we obtain 



and keep only the 0(uj) 



/it? />oc +°° 
^- J dTdT' Tr [T a (E) 



71 — — OQ 



a (2Nt) 2 

g r (E, T)T(E)G a (E, T')] exp(-mu;(T - T'))d E f(E) (20) 
Note that 

/+oo 
dT ^ nuj G r (E,T)exp(-inuj(T -T')) 
n— — oo 

= -2iNr d T >G r (E,T') (21) 
where we have used Eq.(|J). Eq.(po|) now becomes, 

i a = — far [ ^Tv[T a d T g r (T)rg a (T)}d E f , (22) 

T Jo J 2vr 

this is exactly the Brouwer's formulaEEi for adiabatic 
pumping. 

In the following we calculate pumped current using 
Eq.(^|) for a symmetric double <5-barrier structure given 
by U(x) = VqS(x + a) + Vq8{x — a). For this system the 
Green's function G(x,x') can be calculated exactlytH. In 
the adiabatic theory, two pumping potentials are needed 
in order to give a nonzero pumping current. At finite 
frequency, a single pumping-potential was reported to be 
enough to pump a currentlill which is peculiar, and our 
exact non-perturbative theory derived above allows us 
to investigate this situation clearly and unambiguously. 
We chose the single pumping potential to be sinusoidal 
V(x, t) — V p cos(u>t)5(x — xq), and we calculate pumped 
current from the left lead at zero temperature and set 
Vq = 200 unless specified otherwise. A gate voltage v g is 
applied to the double barrier structure so that the reso- 
nant single particle energy level is controlled by it. We fix 
the Fermi level of leads in line with the resonant level at 
v g = 0. Finally the unit is set by h — 2m = q = 2a = 1. 
For GaAs material with well width a = lOOOA, the energy 
unit is E = 0.056mey which corresponds to u> = 13.2 
GHz. The unit for pumped current is 5 x 10~ 10 A. 

Fig.l shows the pumped current versus the gate volt- 
age in the strong pumping regime for xq — —a, V p = 160, 
and uj = 1. Following observations are in order. (1) The 
pumped current displays a series of remarkable plateaus. 
The width of a plateau is equal to hu> whereas the height 
of the plateau depends on V p and Vq in a nonlinear fash- 
ion. (2) The pumped current reverses its direction when 
the resonant level across the Fermi level of the lead, i.e., 
when v g changes sign. As a result, we see from Fig.l that 
the pumped current is quenched near v g — 0. (3) The 



height of the current plateaus for negative v g is consid- 
erablely larger than those of positive v g . To understand 
these results, we rewrite Eq.(@) as 



h = q 



^F n (E)(f n (E) 



f{E)) 



(23) 



with 
F n (E) 



T TT[T L (E)G r (E,E + ntu)r n (E) 



(2Nt) 2 

G a {E + nu, E) - T Ln {E)G r {E + nu, E) 
T(E)G a (E 1 E + nu)] (24) 

We notice that kernel F n consists of two terms: one due 
to photon absorption procesaU and the other to photon 
emission process. These two contributions differ by a 
sign indicating the competition between the photon ab- 
sorption and emission processes similar to that discussed 
in Ref. |ll|. In Fig.2 we plot the integrand of Eq.(||), 
F n (E), for n — 1,2,3 (solid line, dotted line, dashed 
line, respectively). The sidebands due to photon absorp- 
tion (when v g < 0) and photon emission are clearly seen 
(for an n-photon emission process, the energy has been 
shifted by nto due to the transformation from Eq.(^|) to 
Eq.(^3|)). To obtain pumped current, one integrates these 
sidebands over energy. For the n-th sideband (n > 0), the 
integral range is [v g — nu, v g \: it is therefore clear that 
the plateau structure is a direct result of multi-photon 
assisted processes. The positive current is due to pho- 
ton absorption processes and the current reverses its sign 
if photon emission process dominates. The behavior of 
current reversal has lieen observed before for pumping 
in Carbon nanotubescfl, charge quantizationo, and heat 
current generated during the pumpingc3, although from 
completely different origins. Due to the energy depen- 
dence of coupling between the pump and the leads, the 
sideband is asymmetric with a larger peak for photon 
absorption process. As a result, the height of the cur- 
rent plateau is larger for photon absorption process. Our 
numerical plots show that this plateau structure persists 
for different frequencies and pumping amplitudes. Our 
analysis show that there are two effects which affect the 
current plateau structure. The first is barrier heights 
of the quantum well: current plateaus can only be ob- 
served in the strong tunneling regime and they disap- 
pear for low barrier height. The second is temperature: 
the plateau structure is rounded off at finite temperature 
and destroyed when temperature reaches ~ 200m.fr for 
w = 10GHz. 

When the pumping potential is inside of the double- 
barrier structure, we found a similar plateau structure 
due to the same photon assisted processes. In Fig. 3, the 
pumped current as a function of pumping position is 
plotted at a fixed v g = —0.67 for a much smaller pumping 
amplitude v p = 6. We observe that the pumped current 
I p is antisymmetric about the center of the pump. As ex- 
pected for a symmetric structure, I p vanishes identically 



3 



at x = 0. At x = -0.5, we found I p = 2 x 1CT 6 . Fig.3 
shows that I p can increase by three orders of magnitude 
by varying the position of the pumping potential away 
from the barriers. 

In summary, we have derived a non-perturbative and 
exact current expression for parametric pumping at finite 
frequency and finite pumping amplitude. This theory can 
also account for the case of pumping potential located 
in the interior of the scattering region. Our theory re- 
produces all previously known results in both adiabatic 
regime and in perturbative theory. For a double-barrier 
quantum well pump, we predicted current plateaus to 
appear due to multi-photon assisted processes with the 
plateau width given by the pumping frequency. As the 
gate voltage is varied and the resonance level sweeps 
through the Fermi energy of leads, the pumped current 
is found to reverse its direction. 
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FIG. 1. The pumped current versus the gate voltage for 
xo = —a, u = 1, and v p = 160. 

FIG. 2. The integrand of Eq.(23) versus gate voltage for 
n = 1 (solid line) , n — 2 (dotted line) , and n — 3 (dash-dotted 
line). The system parameters are the same as that of Fig.l. 
For illustrating purpose, we have offset the curve by 0.02. 

FIG. 3. The pumped current versus position of the pump- 
ing potential. Here v g = —0.67, u> = 1, and v p — 6. 
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